Soft singularity crossing and transformation of matter properties 
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We investigate particular cosmological models, based either on tachyon fields or on perfect fluids, 
for which soft future singularities arise in a natural way. Our main result is the description of a 
smooth crossing of the soft singularity in models with an anti-Chaplygin gas or with a particular 
tachyon field in the presence of dust. Such a crossing is made possible by certain transformations 
of matter properties. Some of these cosmological evolutions involving tachyons are compatible with 
SNIa data. We compute numerically their dynamics involving a first soft singularity crossing, a 
turning point and a second soft singulatity crossing during recoUapse, ending in a Big Crunch 
singularity. 
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I. INTRODUCTION 

In a recent paper [V\ we have investigated the pos- 
sibility of soft singularity crossing in a model where a 
flat Friedmann universe was filled with dust and anti- 
Chaplygin gas. A soft singularity appears in an expand- 
ing universe when the pressure of the anti-Chaplygin gas 
diverges, causing a — ?► — cx) (a is the scale factor and the 
dot denotes the derivative with respect to cosmic time), 
while a and a remain finite. The energy density of the 
anti-Chaplygin gas vanishes at the singularity while the 
energy density of dust remains finite there, thus implying 
d > 0. The geodesic equations remain regular at the sin- 
gularity and, therefore, they can be continued through. 
This allows the universe to cross the soft singularity. 
Then, a smooth evolution of the universe would require 
further expansion. However, in this case the energy den- 
sity of the anti-Chaplygin gas would become imaginary 
and hence ill-defined. This contradiction and the fact 
that the geodesies can be continued through leads to a 
paradox. In [l| we have solved this paradox by relaxing 
the smoothness condition, leading to the redefinition of 
cosmological quantities as distributions. With this redef- 
inition, it turns out that the universe can revert abruptly 
from expansion to contraction. 

In the present work we study an alternative possibility 
for the continuation of geodesies across the singularity, 
requiring the continuity of the spacetime evolution at the 
expense of certain transformations of matter properties. 

The paper is organized as follows. In Sec. II we dis- 
cuss the sudden singularities which arise in a class of 
fiat Friedmann models such as those describing universes 
filled with the anti-Chaplygin gas without or with dust, 
or driven by a specific tachyon field [2], again without 
and with dust. In Sec. Ill we describe the crossing of 
a soft singularity with accompanying transformations of 
matter in the above mentioned models. Sec. IV presents 



numerical results of the evolution of universes filled with 
a tachyon field and dust, compatible with SNIa data. 
Concluding remarks are presented in Sec. V. We choose 
units c = 1 and SttG/S = 1. 



II. SUDDEN SINGULARITIES IN FLAT 
FRIEDMANN UNIVERSES 

The line element squared of a fiat Friedmann universe 
can be written as 



ds'' = dt^-a^{t)J2{dx''f, 



(1) 



where x" {a = 1,2,3) are spatial Cartesian coordinates. 
The evolution of the universe is governed by the Ray- 
chaudhuri (second Friedmann) equation 



and by the continuity equation for the fluid 
p + 3H{p + p)^0, 



(2) 



(3) 



where, as usual, p is the energy density, p is the pressure 
of matter and H = d/a is the Hubble parameter. The 
first Friedmann equation is 



H' 



P- 



(4) 



Sudden singularities are characterized by finite Hs and 
Hs = — oo (finite dg and dg = — oo) at some finite scale 
factor as- Here, the subscript S denotes the respective 
quantities evaluated at the singularity. These conditions 
can be formulated in terms of energy density and pressure 
of the fluid. The Friedmann ^ and Raychaudhuri ([2]) 
equations show that the total energy density ps is non- 
negative and finite while the pressure diverges p5 = oo . 



It was shown in [3] , 4] and 1] that the geodesies can be 
continued across such sudden singularities as the geodesic 
equations are regular there. The singularity is weak (soft) 
according to the definitions of both Tipler p\ and Krolak 
Q. Although the tidal forces become infinite, the ex- 
tended objects are not necessarily crushed when reaching 
the singularity. 



A. Big Brake 

A special case of sudden singularity is the Big Brake 
singularity, occurring when the energy density vanishes 
at the singularity, ps = ^. 



1. Anti-Chaplygin gas 

One of the simplest models where the Big Brake sin- 
gularity arises is the anti-Chaplygin gas [^. This is a 
perfect fiuid with the equation of state 



P 



A 
P 



(5) 



where A > 0, as opposed to the Chaplygin gas 0, IM| 
which has the equation of state p — —A/ p. The equation 
of state (l5]) arises, for example, in the theory of wiggly 
strings [9]. 

Applied to the anti-Chaplygin gas, the continuity equa- 
tion ^ gives the following dependence of the energy den- 
sity on the scale factor: 



PACh 



B 



-A, 



(6) 



where _B is a positive constant, which determines the 
initial condition. When the scale factor approaches the 
value 



B 



(7) 



during the expansion of the universe, the energy density 
of the anti-Chaplygin gas vanishes, and its pressure grows 
to infinity. Accordingly, the deceleration also becomes 
infinite. 

As it was shown in [l| , after crossing of this singularity 
the universe starts contracting towards a Big Crunch. 



2. The tachyon field with trigonometric potential and 
transition to a Bom-Infeld type pseudo-tachyon field 

A Big Brake singularity was first found in a specific 
tachyon model introduced in ^ . The Lagrangian density 
of a tachyon field is [10| 



L 



-V{T)^l-g^Hd,T){d,T), (8) 



where V{T) is a potential. For a spatially homogeneous 
field T{t) the expression ([8|) becomes 



L^~V{T)Jl-gt*T^ 



(9) 



This field corresponds to an ideal fiuid with energy den- 
sity 



Pt 



V{T) 



and pressure 



Pt 



^-V{T)Vl-T^. 



(10) 



(11) 



The Lagrangian density as well as px and px are well 
defined for T^ < 1. The field equation is 



T 



3HT 



Vi 



= 0. 



1-T2 ■ """ ' V 
The following potential was studied in 0: 



with 



ViT) = 



y — cos 



Av/l-(l + %2 

1-y^ 



-v/A(i + fc)r 



(12) 



(13) 



(14) 



where A > and — 1 < A: < 1 are model parameters. The 
case fc > is of particular interest, because it reveals two 
unusual features: a self-transformation of the tachyon 
into a pseudo-tachyon field and the appearance of a Big 
Brake cosmological singularity. For fc > 0, the potential 
([T ^ - P^ . is well-defined in the range 



where 



y* <y <y*, or r4 > T > T3, 

y, = (l + fc)-l/2^ 



T. 



3^A(l + fc) 



arccos(l + k) 



^1/2 



Ta 



i^K[i + k) 



7r-(l + /e)-i/2 



(15) 

(16) 
(17) 

(18) 



Note that the dynamical system is invariant under the 
simultaneous change 



y 



-y, T 



^T 



(19) 



Since V^ > the pressure is negative allowing for an ac- 
celerated expansion of the universe. When they reach 
the attractive critical point (y = 0,T = 0) the trajec- 
tories correspond to an exact de Sitter expansion of the 



universe. 



The lines T = ±1 (with the exception of the 
corner points (±y,, T = ±1)) in the {T,T) space corre- 
spond to a standard Big Bang singularity (see Fig. 5 of 
[J|, which reproduces Fig. 4 in 0]). 

However, some trajectories can reach the corner points 
where the geometry is not singular. Hence, the trajec- 
tories can be continued across these corner points, be- 
yond which |Tp becomes larger than 1. The potential V 
and the kinetic term in the Lagrangian density © be- 
come imaginary accross the corner points, however their 
product remains real. Thus for T^ > 1 the correct La- 
grangian density (describing a Born-Infcld type pseudo- 
tachyon field) is 



L 



W{T)sfs 



*T2 - 1, 



where 



W{T) = 



Av/(l + fc)y2-l 



l-y2 

This Lagrangian is well-defined in the ranges 
- 1 < y < -?/*, or Tmax >T >Ti 
and 



(20) 



(21) 



(22) 



y, < y < 1, or Tg > T > 0, 



with 



T,„ax = 7r/3v/A(my. 
The energy density and pressure are now 

W{T) 



Pt 



Vt^ - 1 ' 



Pt 



w{t)Vt^ 



1. 



(23) 



(24) 



(25) 



Since the pressure is positive, the expansion of the uni- 
verse is slowing down. The field equation is 



T 



1-T2 



5HT 






0. 



(26) 



The universe runs into a soft singularity somewhere in 
the ranges ^<j/<lasT— >■ — cjo or — 1 < y < — y» 
as T -J> oo 0. From Eqs. (I1I|)-(I1S|), the potential W 
is finite and pT ^ 0, Pt — >■ oo at the soft singularity. 
Equivalently, this means that Hg — while Hs = — oo. 
It was shown in [4| that the evolution of the universe can 
be continued accross the singularity, the universe starts 
recoUapsing and eventually ends in a Big Crunch singu- 
larity. 



B. Introducing a dust component 

The dust is a perfect fluid with vanishing pressure, 
whose energy density is 



Pm,0 



(27) 



where Pm,o is a positive constant, characterizing the 
quantity of matter in the universe today (ag = 1). There- 
fore, if a cosmological model with dust evolves into a 
sudden singularity the energy density of dust-like matter 
remains finite ipm)g > 0. Then, the Hubble parameter 
does not vanish at the singularity as in the case of the 
Big Brake. This makes it more difficult and delicate to 
describe what happens after reaching a soft singularity. 



1. Anti-Chaplygin gas 

A soft singularity arising in a two-fluid model contain- 
ing dust and anti-Chaplygin gas was investigated in [l|. 
The Hubble parameter is positive at the singularity, re- 
quiring a further expansion of the universe. Then a para- 
dox arises: if the universe continues to expand, beyond 
the singularity the expression under the sign of the square 
root in Eq. ([5]) becomes negative and the energy density 
of the anti-Chaplygin gas becomes ill-defined. 

A mathematically consistent way out of this situation 
is an abrupt replacement of the cosmological expansion 
by a contraction at the price of introducing distributional 
cosmological quantities [1] . 

In the next section, we investigate an alternative pos- 
sibility requiring the smoothness in the evolution of the 
Hubble parameter but allowing for a change in the equa- 
tion of state ([5]) . 



2. Born-Infeld type pseudo-tachyon field with trigonometric 
potential 

In the model suggested in Q] the Born-Infeld type 
pseudo-tachyon field runs into a soft Big Brake singular- 
ity at some point during the expansion of the universe. 
What happens however in the presence of a dust compo- 
nent? Does the universe still run into a soft singularity? 

In order to answer this question, we rewrite Eq. (j26p 
as 



t= (T^- 1) ( 3iJr- 



W.T 



(28) 



In the left lower and in the right upper stripes (see Fig. 
4 of [2]), where the trajectories describe the expansion 
of the universe after the transformation of the tachyon 
into the pseudo-tachyon field, the signs of T, of T and of 
the term W^t/W coincide. A detailed analysis based on 
this fact was carried out in ^ and led to the conclusion 
that the universe encounters the singularity as T — > Tg 



{Ts > or Ts > Tmax) , |T| ^- 00. The presence of dust 
cannot alter this because it increases the influence of the 
term 3HT, and hence, accelerates the encounter with the 
singularity. Indeed, consider two trajectories, crossing 
one of the corners (i.e., undergoing the tachyon-pseudo- 
tachyon transition) under the same angle in phase space 
(cf. Fig. 4 in [2]), one in the absence of dust, the other in 
the presence of dust. For both trajectories the signs of T 
and of T in (|28l) coincide and the increase that the value 
of H undergoes when dust is present makes the growing 
of \T\ more abrupt. On the other hand, the evolution of 
the tachyon field, approaching the corner point is slowed 
down by the presence of dust, because, in this case (inside 
the rectangle of the phase space) (T^ — 1) is negative 
and therefore T and T have opposite signs. Summing 
up, we may say that the presence of dust accelerates the 
evolution of the pseudo-tachyon whereas it slows down 
the evolution of tachyon. 

What is important is that the presence of dust changes 
in an essential way the time dependence of the pseudo- 
tachyon field close to the singularity. Indeed, as it was 
shown in Q , in the absence of dust one has 



T = Tbb± 



( 4 
\3W{Tbb) 



1/3 



[tBB ~ t) 



1/3 



(29) 



(see Eq. (29) in jj]). The upper (lower) sign corresponds 
to the left lower (right upper) strip in Fig. 4 of [2], where 

Ymit^tBB ^ — ^°^ ( linit^tBB T — oo). In the presence 
of dust one has, instead, 



T^Ts± 



3Hs 



Vt^^t, 



(30) 



where Hs is the nonvanishing value of the Hubble pa- 
rameter given by 



Hs 



Pm,0 



(31) 



Here we have taken advantage of the fact that in Eq. 
the terms 1 and W^t/W can be neglected with respect to 
T^ and 3HT, respectively. It is easy to see that a smooth 
continuation of expression (I30|) is impossible in contrast 
to the situation without dust (f^ . 

Thus, the presence of dust is responsible for the ap- 
pearance of similar paradoxes in both the anti-Chaplygin 
gas and tachyon models. 



III. CROSSING THE SOFT SINGULARITY 
AND TRANSFORMATIONS OF MATTER 

As mentioned earlier (see Introduction of the present 
paper and the Concluding remarks in [l|) the mathemat- 
ically self-consistent scenario, based on the treatment of 
physical quantities as generalized functions and on the 
abrupt change of the expansion into a contraction, may 



look counter-intuitive from the physical point of view. 
Indeed, such a behavior displays features which are anal- 
ogous to the phenomenon of the absolutely elastic bounce 
of a hard ball from a rigid wall, as studied in classical me- 
chanics. In the latter case, it is the velocity and the mo- 
mentum of the ball which change their direction abruptly. 
Hence, an infinite force acts from the wall onto the ball 
during an infinitely small interval of time. 

In reality, the absolutely elastic bounce is an idealiza- 
tion of a process taking place in a finite, though small, 
time-span, during which inelastic deformations of the ball 
and of the wall occur. This implies a more complex and 
realistic description of the dynamical process of inter- 
action between the ball and the wall. Hence, we are 
naturally led to assume that something similar should 
occur also in the models of an anti-Chaplygin gas or 
a tachyon whenever dust is present. We expect that 
the smoothing of the process of the transition from an 
expanding to a contracting phase should include some 
(temporary) geometrically implied change of the equa- 
tion of state of matter or of the form of the Lagrangian. 
We know that such changes have been considered in cos- 
mology. For example, a tachyon-pseudo-tachhyon trans- 
formation, driven by the continuity of the cosmological 
evolution, took place in the tachyon model [2] (see also 
subsection II. B. 2 of the present paper). In a cosmologi- 
cal model with the phantom field with a cusped poten- 
tial [12l | , transformations between phantom and standard 
scalar field were considered. Thus, it is quite natural to 
assume that the process of crossing of the soft singularity 
should imply similar transformations. 

However, the situation is now more complicated. It is 
not enough to require the continuity of the evolution of 
the cosmological radius and of the Hubble parameter. It 
is also necessary to make some hypotheses about chang- 
ing the equation of state of matter or the form of the 
Lagrangian. 

We solve the problem as follows. Considering first 
the anti-Chaplygin gas with dust, we require a minimal 
change in the form of the dependence of the energy den- 
sity and of the pressure on the cosmological radius, upon 
crossing the soft singularity. This will require replace- 
ment of the anti-Chaplygin gas with a Chaplygin gas with 
negative energy density^. Next, we consider the cosmo- 
logical model based on a pseudo-tachyon field with con- 
stant potential and in the presence of dust. It is known 
that the energy-momentum tensor for such a pseudo- 
tachyon field coincides with that of the anti-Chaplygin 
gas (relating the Chaplygin gas to the tachyon field with 
constant potential was considered in [iJI). We derive how 
the pseudo-tachyon Lagrangian transforms using its kin- 
ship with the anti-Chaplygin gas. In this way, we arrive 
at a new type of Lagrangian, belonging to the "Born- 



^ A Chaplygin gas with negative energy density has been consid- 
ered earlier ^13| in a different context. 



Infeld family" . Finally, we extend this transformation to 
the case of the trigonometric potential. 



A. Anti-Chaplygin gas 

It follows from Eqs. ([5]) and ([6|) that the pressure of 
the anti-Chaplygin gas 



A 



(32) 



tends to +(X) when the universe approaches the soft sin- 
gularity, e.g. when the cosmological radius a ^ as (see 
Eq. ([7])). Requiring the expansion to continue into the 
region a > as, while changing minimally the equation of 
state, we assume 



A 



P 



or, in other words, 



A 



P = 



A-4 



^1 



for a > as- 



(33) 



(34) 



Thus, the pressure passes through -|-c» conserving its 
sign, thus providing in such a way the continuity of the 
cosmological evolution. It is crucial that p does not 
change sign in order to keep a decelerated expansion. 
The energy density p evolves continuously, and so does 
its derivative with respect to volume. Combining (p4)) 
with the energy conservation law ([3]) we obtain 



p=-\lA -, for a > as. 



(35) 



so that for a > as the energy density and the pressure 
satisfy the Chaplygin gas equation of state 



P 



A 

P' 



(36) 



Therefore, at the singularity crossing, the anti-Chaplygin 
gas transforms into a Chaplygin gas with negative energy 
density. After crossing of the singularity the Friedmann 
equation is 



H2 = PI!1^-VA^Ii-[^) 



(37) 



and it follows from Eq. ((37|) that, after achieving the 
point of maximal expansion a = flmax, where 



Pmfl 

^4~ 



1/6 



(38) 



the universe begins contracting. During this phase, as 
it achieves again a = ag, it stumbles once more upon a 



soft singularity, whereupon the Chaplygin gas transforms 
itself back into anti-Chaplygin with positive energy den- 
sity and the contraction continues until hitting the Big 
Crunch singularity. 

Whereas in [1] we envisaged an abrupt change from 
expansion to contraction through the singularity, with 
a jump in the Hubble parameter, we show here that a 
continuous transition to the collapsing phase is possible 
if the equation of state of the anti-Chaplygin gas has 
some kind of "phase transition" at the singularity. 



B. Pseudo-tachyon field with a constant potential 

For a pseudo-tachyon field with constant potential 
W{T) = Wo, the energy density (p4l) and the pressure 
(PSj) satisfy the anti-Chaplygin gas equation of state ^ 
with 



A = Wl 



(39) 



Solving the equation of motion for the pseudo-tachyon 
field (HH) with W{T) = Wo, one finds 



rp2 _ 1 



1 



(4) 



(40) 



and we see that a soft singularity arises at a = a^ with 
T^ -^ +0O. 

The new Lagrangian, which gives the correct energy 
density and pressure satisfying a Chaplygin gas equation 
with negative energy density is 



L = Wo^5«T2 + l, 



a> as 



givmg 



and 



P = WoVt^Ti 



W. 



P = 



VT^+l 



(41) 



(42) 



(43) 



Lagrangian (j41l) characterizes a new type of Born-Infeld 
field, which we may call "quasi-tachyon" . 

For an arbitrary potential the Lagrangian reads 



with equation of motion 



a> as (44) 



T ■ Wt 



T2 -t- 1 ■ """ W 

and energy density and pressure are respectively 

W{T) 



P = 



Vf^ + i 



(45) 



(46) 



and 



p = W{T)\/t^ + l. 
If W{T) = Wo, the solution of equation gS]) is 

1 



rp2 



(4) 

and the energy density evolves as 



6 ' 

-1 



(47) 



(48) 



Pt 



*""-(t 



(49) 



The evolution of the universe coincides with that of a 
universe with anti-Chaplygin gas and dust. 

Thus, the transformation from anti-Chaplygin to 
Chaplygin gas with negative energy density corresponds 
to a transition from a pseudo-tachyon field with La- 
grangian (^0)) with constant potential W{T) = Wq to 
a new Born-Infcld type quasi-tachyon field, with La- 
grangian (l4Tt . 



The tachyon model with trigonometric potential 
and dust 



In the vicinity of the soft singularity, it is the "fric- 
tion" term 3HT in the equation of motion (pS)) . which 
dominates over the potential term W^t/W. Hence, the 
dependence of W{T) on its argument is not essential and 
a pseudo-tachyon field approaching this singularity be- 
haves like one with a constant potential. Thus, it is rea- 
sonable to assume that upon crossing the soft singularity 
the pseudo-tachyon transforms itself into a quasi-tachyon 
with Lagrangian (H^ for any potential W{T). 

We now study the dynamics of the model with trigono- 
metric potential p^ - (ITi|) in the presence of dust. 

The behavior of the quasi-tachyon field close to the 
soft singularity can be derived from Eq. (j45p in the same 
way as the corresponding behavior of the pseudo-tachyon 
field derives from Eq. (pS)) . In analogy with Eq. (150]) . 
we obtain the quasi-tachyon behavior 



T^TsT 



3Hs 



Vt^^ 



(50) 



and the two formulas match with each other through the 
singularity. 

In order to analyze the dynamics of the field in the 
presence of dust, it is convenient to concentrate ourselves 
on the processes as they occur, say, in the left lower strip 
of the phase diagram of the model, to facilitate compar- 
ison with earlier studies of the tachyon model dynamics 
without dust in @, S] ■ The relative signs in the equations 
of motion of the term with the second derivative T and of 
the friction term 3HT are opposite for pseudo-tachyons 
and quasi-tachyons. This means that after crossing the 



soft singularity the time derivative T grows while its ab- 
solute value decreases. At the same time the value of T 
is decreasing while the potential W{T), given by (|2T|) is 
growing. 

Hence the absolute value of the negative contribution 
to the energy density of the universe induced by the 
quasi-tachyon grows while the energy density of the dust 
decreases due to the expansion of the universe. Thus, 
at some moment the total energy density vanishes and 
the universe reaches the point of maximal expansion, af- 
ter which the expansion is replaced by a contraction and 
the Hubble variable changes sign. The change of sign 
of the friction term 3HT implies the value of T to de- 
crease and at some finite moment of time the universe hits 
again the soft singularity when T — >■ — oo. Upon cross- 
ing this singularity the quasi-tachyon transforms back to 
pseudo-tachyon and the relative signs of the terms with 
the second and first time derivatives in the equation of 
motion change once again. After this, the time deriva- 
tive of the pseudo-tachyon field begins to grow and the 
universe continues its contraction until it hits the Big 
Crunch singularity. 

It was shown in Q that, for the case of the tachyon 
model with trigonometric potential and without dust, the 
encounter of the universe with the Big Crunch singularity 

occurs at T = and T = —J i-^ . One can show that the 

presence of dust does not change these values. Indeed, 
consider the behavior of the pseudo-tachyon field when 

T — > 0, T — > ^y-^^- It follows from the expressions 

p4| and (|25)) that the ratio between pressure and energy 
density behaves as 



p k 



(51) 



i.e. in the vicinity of the Big Crunch singularity the 
pseudo-tachyon field behaves as a barotropic fluid with 
the equation of state parameter t: > 1. This means that 
the energy density of the pseudo-tachyon field grows as 



P' 



1 



a3(i+i) 



(52) 



with a -T> 0, namely much more rapidly than the dust 
energy density. Thus, one can neglect the contribution 
of dust in the regime of approach to the Big Crunch sin- 
gularity and the description of the evolution of the uni- 
verse to this point coincides with that of the pure tachyon 
model 0. 



IV. FUTURE EVOLUTION OF THE TACHYON 

FIELD WITH TRIGONOMETRIC POTENTIAL 

AND DUST: NUMERICAL RESULTS 

The tachyon model with trigonometric potential was 
tested in [ll| by comparing it with SNIa data. In that 
paper we found the range of values of the model param- 
eter k and tachyon field initial conditions fitting well the 



SNIa data. Then we studied future evolutions starting 
from acceptable initial conditions. While a subset of the 
corresponding trajectories leads to a de Sitter expansion, 
a complementary subset of trajectories leads to a Big 
Brake singularity. The evolution after the Big Brake sin- 
gularity crossing was described in [J]. 

In subsection IVA we investigate the compatibility 
with SNIa data of this dark energy model in the pres- 
ence of dustlike matter. We use the Union2 SNIa data 
set [ial ■ We show that the model fits the SNIa data well 
also in the presence of dust like matter. 

In subsection IVB we investigate the future evolution 
numerically for those trajectories which run into the soft 
singularity at time ts^- We give specifically the time 
intervals measured from today for the following events: 
tdec when the cosmic expansion becomes decelerated; t* 
corresponding to the tachyon-pseudo-tachyon transfor- 
mation (crossing the corner of the rectangle in the phase 
portrait of the model); tg-^ when the first soft singularity 
is reached; ttum corresponding to the turning point when 
the universe starts contracting; i^^ when the second soft 
singularity is reached; and finally the time tsc of the Big 
Crunch. 



through the values {0.03, 0.09, 0.15, 0.21, 0.27, 0.33}. 
As in paper [11[ we avoid the double coverage of the 
parameter space (the model has a symmetry given by 
Eq. (IT51) ) by replacing Tq 11| with the new variable'^: 



1 



Xo 



rp2 

^0 



(55) 



The initial conditions xo and j/q vary inside the rectangle 
^ < Xo < l,|yo| < 0.845. Finally, we introduce the 
luminosity distance cIl whose evolution is given by 



d_ 

dz 



1 + 2 



1 



(56) 



Fitting to the supernovae data involves a x test, as 
described in Refs. [ia|, [ll|. In Fig [T] we show the y^ 
values in the parameter plane of the initial conditions 
(?/0i xq). The contours correspond to the 1 a, respectively 
2cr, confidence levels with x^ = 570.34 and respectively 
X^ = 612.33. 



B. Future evolution 



A. Test with supernovae data 

The tachyon field violates the strong energy condi- 
tion when T^ < 1, as required by a dark energy can- 
didate. For a reasonable fit with supernova data we as- 
sume T^ < 1 . In the regime where a varies monotonically 
with time it may be convenient to replace the cosmolog- 
ical time with a monotonic function of the scale factor 
as a new independent variable. We choose the redshift 
1 -\- z = ao/a as a new independent variable^ (here and 
henceforth the subscript refers to the value of the re- 
spective quantities at the present epoch). 

The model depends on the parameters k and A. For 
given values of these parameters the possible solutions 
depend on the quantity of dust and on the initial con- 
ditions hojTq for the tachyon field. However, the Fried- 
mann equation implies the following constraint 



To 



ni,[i-ii + k)yi] 



(1 - a„,o)' (1 - yg) 

where (remembering our convention 87rG/3 — 1) 

Pm,0 



A 



f^A,0 — TTfj f^m,0 



Hi' 



(53) 



(54) 



showing that A is determined by the values of the other 
parameters. In what follows we fix fc = 0.4 and vary Qm,o 



As in the preceding papers [j,lll| we study numerically 
the future evolution of the universe starting with initial 
conditions compatible with SNIa data. However, our task 
now is technically more complicated due to the presence 
of dust. As a matter of fact, we shall have to consider five 
different regimes, where different systems of dynamical 
equations are used and we should provide four accurate 
matching between these evolutions. First, the universe 
starts its evolution at somepoint in the rectangle on the 
phase space of Fig. 4 of I3|. Here the field T satisfies 



the equation of motion p2p and the right-hand side of 
the first Friedmann equation includes the contribution of 
dust ((27l) and of the tachyon field ([10]) . After the crossing 
of the corner (at t*), the tachyon field transforms into 
a pseudo-tachyon field with equation of motion (pS)) and 
energv densitv (|24p . This is the second regime. The third 
regime enters into action after the first crossing of the soft 
singularity (at is^), when the pseudo-tachyon transforms 
itself into a quasi-tachyon with equation of motion ps)) 
and energy density (HS)) . After the passing of the point of 
maximal expansion of the universe (at ttum) we enter into 
the fourth regime when the universe starts contracting. 
After the second soft singularity crossing (at ts^ ) we have 
the fifth regime, where the quasi-tachyon converts itself 
again into a pseudo-tachyon. Finally, the universe ends 
in a Big Crunch (at tsc)- The corresponding times are 
shown in Table [J for 0^,0 — 0.03 and in Tabic |TT] for 
^m.o = 0.27. 



^ We note that since the Friedmann equation is a first integral it 
can be used as a check of the accuracy of the numerical integra- 
tion. 



^ The parameter xq is denoted by too in [ll| . 
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FIG. 1; (Color online) The fit of the luminosity distance vs. redshift for the parameters k — 0.4 and flm = 0.03 (upper left), 
0.09 (upper right), 0.15 (middle left), 0.21 (middle right), 0.27 (lower left), 0.33 (lower right), in the parameter plane (j/o, xo) 
in the range \yo\ < 0.845 where the potential V is well defined. The contours refer to the 68.3% (Icr) and 95.4% {2a) confidence 
levels. The colour code for x. is indicated on the vertical stripes. The clear tendency with increasing the dust component is 
that the parameter xo approaches its maximally allowed value (representing To = 0). Higher values of i},n render the fit with 
the supernovae outside the la region. 



r 



These times have been computed assuming for Hq the 
value 70 km s~^ Mps ~^. It is known that there is a cer- 
tain discrepancy between the value of the Hubble param- 
eter arising indirectly from the cosmic microwave back- 
ground and baryon acoustic oscillations [13] , and the one 
more directly obtained from local measurements of the 



relation between redshifts and distances to sources [l8| 
(for a recent analysis of this problem see [.19|). The for- 
mer gives H^^^ = 67.89 ± 0.77 km s^^ Mps^^, while 
the latter gives iJ^°™' = 73.8 ± 2.4 km s^^ Mps^^. Nev- 
ertheless, the precise value of Hq is not so important for 
our study, hence, we have taken an intermediate value. 



TABLE I: Key times in the evolution of tachyon universes for k — 0.4 and Qm.o = 0.03 are given. The first two columns give 
initial values xq, yo in agreement with supernovae data at the la confidence level. Remaining columns starting from left give 
the successive times (measured from the present time): tdec when the expansion becomes decelerated, the tachyonic tansition 
time f«, the first soft singularity crossing time ts-^, the turning point ttum, the second soft singularity crossing time ts2 and 
finally the Big Crunch time tsc- Times are given in 10 yrs unit and calculated assuming Ho = 70 km s^^ Mpc~ . 
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TABLE 11: Same as in Table H for H^.o = 0.27. 
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tSi 


ttum 
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0.65390 
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7.81302 


7.81738 


8.15389 


-0.65 


0.995 


21.07994 


21.33077 


22.22677 


22.22681 


22.22689 


22.57314 



Now we can turn to the analysis of the Tables I and II. 
In Table U different times measured from today are given 
for a low amount of dust ^m,o — 0.03 and in Table HIl for 
^m,o — 0.27. So we see the effect of the addition of dust 
in a systematic way. Three comments are in order here. 

First, the time interval i* between today and the first 
transition into the pseudo-tachyon varies considerably 
within the set of trajectories compatible with the su- 
pernovae data. Namely, it varies from 0.9 to 9.5 billion 
years for Q.mfl — 0.03 and from 1 to 21 billion years 
^m,o = 0.27. Second, the time intervals between t.^, and 
the Big Crunch time Ibc are practically constant (about 
1.1 billion years for the first case and about 1.2 billion 
years for the second case) . A similar property was found 
in the model without dust [4|. Third, the time interval 
between the two soft singularity crossings 1^2 — tsi de- 
creases strongly (from 810 thousand years to 70 thousand 
years for the first case and from 0.03 to 0.0002 billion 
years for the second case) when the value of t, increases. 
This can be ascribed to the fact that the density of dust, 
at the moment of the first soft singularity crossing ^5^ , 
for the universes with high values of t* is greatly reduced 
compared to those with small values of i* . Indeed, in the 
absence of dust the two values tsi and tg^ coincide and 
we have a unique Big Brake singularity. 

On Figure 2 the evolutions are shown in the three- 
dimensional coordinate space x,y,z for six different val- 
ues of ftm,o- For the trajectories ending in a de Sitter 
space, the final point has coordinates (1, 0,-1). For other 
trajectories we present only the evolutions until the first 
soft singularity crossing. Generally, the sets of initial con- 



ditions, compatible with the supernovae data (the regions 
in the plane (x, y) at z = 0) decrease as the quantity of 
dust increases and vanish for ftm,o > 0.33. Also, as ^m,o 
increases, the number of trajectories going to a soft sin- 
gularity is decreasing compared to those ending in a de 
Sitter space. 

This work is done with the spatially-flat paradigm in 
mind. However, as this model is constrained using SNIa 
data, it is interesting to relax the assumption of flatness 
in this case and to consider also non-flat universes with 
a spatial curvature allowed by observations. Indeed in a 
spatially closed universe the curvature and matter terms 
in the Friedmann equation could cancel each other at 
some (negative) redshift z^urv 



l + Zc 



\n 



fe,0| 



a 



m,0 



(57) 



The quantity ^k.o is strongly constrained by observa- 
tions, -0.0065 < flkfi < 0.0012 (95% C.L.) with central 
value rikM = —0.0027 [22|. Hence a slightly spatially 
closed universe is favoured. 

Of course the tachyon, like any scalar field model and 
in sharp contrast to the anti-Chaplygin gas, does not have 
a barotropic equation of state. Therefore the amount of 
expansion needed to reach the soft singularity depends on 
the initial conditions. It is quite clear however that for 
models studied here we will have jzsj < Izcttrt)!- Hence 
the kind of problem considered in this paper, and the 
mechanism suggested in order to cross the soft singular- 
ity, will remain even in the presence of a tiny curvature. 
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FIG. 2: (Color online) The future evolution of those universes, which are in a 68.3% confidence level fit with the supernova 
data. The la contours (black lines in the z — plane) are from Fig [T] (the parameter plane {yo,xo) is the 2 = plane here). 
The sequence of figures and the values of Qm,o are the same as on Fig. [l] The point (1,0, —1) is the de Sitter final state. 



But we conjecture that peculiar initial conditions do exist 
for which this is no longer the case. 



V. CONCLUDING REMARKS 



Soft cosmological singularities known since the eighties 
[20] ■ have been attracting growing attention during the 
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last few years [2i|. In this paper we have continued the 
investigation of particular cosmological models based on 
tachyon fields or perfect fluids (introduced in paper [2]), 
for which soft singularities arise in a natural way. The 
main result of our investigation is the description of a 
smooth crossing of soft singularities, arising in models 
with anti-Chaplygin gas or of a particular tachyon field in 
the presence of dust. Such a crossing is accompanied by 
certain transformations of matter properties, embodied 
in a change either of equation of state or of Lagrangian. 

The interesting feature of the tachyon model is that 
there exist cosmological evolutions whose past is com- 
patible with the supernova data and whose future reveals 
"exotic phase transitions" which are here described in de- 
tail. We have performed a detailed numerical analysis of 
these evolutions. 

All our studies, both theoretical and numerical, were 
performed assuming a a spatially-flat universe. Next in- 
teresting step for the study of dark energy models pos- 
sessing soft future singularities is the inclusion of spa- 
tially closed universes. Indeed, observations do allow 
for a tiny spatial curvature, a positive curvature being 
slightly preferred. While a tiny viable curvature will not 
change the situation for most models studied in this pa- 
per, a larger number of situations can arise in the pres- 
ence of spatial curvature for the tachyon models because 
of their rich dynamics. Indeed, if the universe reaches the 
point of maximal expansion before occurence of the soft 



future singularity, the latter will not occur at all. In the 
case of our tachyon model this can happen for specific 
initial conditions. If for some peculiar initial conditions 
the turning point and the soft singularity coincide the 
latter retains its character of a Big Brake singularity. (In 
another dark energy model, based on a standard scalar 
field, such an interplay between turning point and the 
encounter with a soft singularity was considered in 23] ) . 
For a comprehensive investigation of these situations a 
more detailed study is required, both theoretical and nu- 
merical and this is left for future work [2J]. In contrast, 
the possible situations in the case of the anti-Chaplygin 
gas are more straightforward. 

Another interesting direction of development of the 
present work is the consideration of cosmological pertur- 
bations and their possible influence on the structure of 
sudden singularities and on the conditions of their cross- 
ing. To our knowledge no systematic study of this kind 
appeared yet in the literature. 
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